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Ford Circles and Spheres 


Sam Northshield* 


Abstract 

Ford circles are parameterized by the rational numbers but are also the 
result of an iterative geometric procedure. We review this and introduce 
an apparently new parameterization by solutions of a certain quadratic 
Diophantine equation. We then generalize to Eisenstein and Gaussian 
rationals where the resulting “Ford spheres” are also the result of iterative 
geometric procedures and are also parameterized by solutions of certain 
quadratic Diophantine equations. We generalize still further to imaginary 
quadratic fields. 


1 Introduction 

The set of Ford circles form an arrangement of circles each above but tangent 
to the x-axis at a rational number, with disjoint interiors, that is maximal in 
the sense that no additional such circles can be added. 




Figure 1: Ford circles and their geometric construction 

They provide a natural way of visualizing the Diophantine approximation 
of real numbers by rationals, are used in the “circle method” of Ramanujan 
and Hardy, and give a geometric way of looking at continued fractions. They 
form part of an Apollonian circle packing, a topic of intense current study, as 
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well as provide some connection to some older open problems (e.g., see m 
for connection to the Riemann hypothesis, [6] for Hausdorff dimension). Three 
dimensional analogues are then certainly of interest. 

In Section 2, we construct the set of Ford circles in three different ways. First, 
to each rational number a/b (in lowest terms), assign a circle above but tangent 
to the x-axis at a/b with radiusl /2b 2 . The set of circles thus parameterized by Q 
is the set of Ford circles which we denote by V. A geometric construction of this 
set starts with circles of radius 1/2 above but tangent to the a?-axis at each of 
the integers. We then proceed inductively: given any two circles tangent to each 
other, we add the unique circle between and tangent to those two (see Figure 
1). The maximal collection formed by this procedure is denoted Q. Third, for 
every relatively prime integer solution (a, b, c) of a 2 + b 2 + c 2 = (a + b + c) 2 
we define a new a circle above and tangent to the x-axis at b/(a + b) with 
radius l/2(a + b). In this last construction, the tangent point has “projective 
barycentric coordinates” (a, b) and so we call this last construction barycentric 
and denote the collection by B. These three collections of circles are shown to 
be all the same. 

In Section 4, we define a family of Ford spheres and three parameterizations 
of it that are analogous to the constructions of Ford circles. For w the cube root 
of unity with positive imaginary part, we may parameterize a family of spheres 
by Eisenstein rationals Q(w), as we did for Ford circles: given two relatively 
prime Eisenstein integers a, /3, we construct a sphere in C x R + tangent to C at 
a/(3 with radius 1/2|/3| 2 . A geometric construction of this family of spheres is 
then given. Starting with spheres of radius 1/2 at every point in the triangular 
lattice {m + nijom, n € Z}, iterate the process: given any three mutually tangent 
spheres, add the two uniquely determined spheres each of which is tangent to 
the original three and to the complex plane. This process is “tetrahedral” in the 
sense that for every three mutually tangent spheres (i.e., their contact graph is a 
triangle), there is a fourth so that the contact graph of all four is a tetrahedron. 
Finally, the relatively prime integer solutions of ( a + b+c+d ) 2 = a 2 + b 2 + c 2 +d 2 
parameterize these spheres barycentrically. 

In Section 5, we define another family of Ford spheres and three parame¬ 
terizations of it that are based on the Gaussian integers. That is, we construct 
family of spheres as in Section 3 using relatively prime Gaussian integers in¬ 
stead of Eisenstein integers. There is a corresponding geometric construction, 
this time based on octahedra instead of tetrahedra. Finally, the barycentric pa¬ 
rameterization is in now terms of solutions of (a+b+c+d) 2 = 2 (a 2 +b 2 +c 2 +d 2 ). 

In Section 6, we consider a definition of Ford spheres parameterized by Q(<r) 
where Z[cr] = for some other positive integers D besides 1 and 3 (so 

that Z[er] is a UFD). In this situation, we develop a barycentric construction as 
well. In the three cases where Z[a\ is a Euclidean domain that are not covered by 
Sections 3 and 4, we conjecture a geometric construction for the corresponding 
Ford spheres. 

In Section 7, we present a summary of the paper and some directions of 
further research. 
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The Ford spheres of associated with uj (Section 4) have been considered by 
Rieger [Ml GS], and those associated with i (Section 5) by Pickover [23] ■ There 
is doubtless some overlap with these papers as well as with the material covered 
in a series of papers by Asmus Schmidt (23 E3 [2H E2 • These last four papers 
are concerned with Diophantine approximation and geometry relate to the cases 
where D = 1,3, 2,11 respectively. Sullivan m deals with quite general arrays 
of spheres C x R + tangent to C. 

Some of the novelties that appear in this paper are as follows. The “barycen- 
tric” parameterizations of Sections 2,4,5 and 6 in terms of solutions of Diophan¬ 
tine equations (e.g.,Equations 7, 8 , 11) is new. We show that for such solutions 
( a,b,c,d ) of (7), \a + b + c\ = | 7| 2 for some Eisenstein integer 7 [Cor. 4.6]. Sim¬ 
ilarly, for solutions (a, b, c, d) of ( 8 ), |a + b\ is a sum of two integer squares and 
|a + 6 + c| is the sum of norm squares of two Eisenstein integers [Cor. 5.10]. We 
introduce a quadratic form on pairs of solutions of equation 7 which takes on 
the value 1 if and only if the corresponding spheres are tangent [Th. 3.4]. Also 
in Section 3, we are able to redefine the Poincare extension of a Mobius trans¬ 
formation in terms of its action on spheres [Prop. 3.2]. In Section 6 , the Ford 
spheres parameterized by Q(<r) are introduced. It is shown also these spheres 
can be parameterized barycentrically (in terms of solutions of equation ( 11 )). 
An intriguing connection is made between these solutions and a group related 
to “secant addition” of m- 

2 Ford circles 

We say a circle in the x, y-plane is normal if it is above and tangent to the 
axaxis. For fgR and r > 0, let C(t, r) be the circle with center ( t , r) and radius 
r. Hence, C(t,r) is normal. We note that every normal circle can be uniquely 
represented as C(t,r) for some t,r. By the Pythagorean theorem, two circles 
C(t,r) and C(t',r') are tangent (we write C(t, r)\\C(t', r’)) if and only if 

(t - t') 2 + (r — r') 2 = (r + r') 2 

or, equivalently, 

(t — t’) 2 = 4 rr'. (1) 

Given o,i£l with b > 0, we define 

Then every circle above and tangent to the x-axis can be uniquely represented 
as C a ,b for some real a,b (b > 0 ): 

C(t,r) = C t / y/ ^ 1 i^p. 

By (1), two such circles are tangent, i.e., C at b\\C Ct d, if and only if |ad — 5c| = 1. 
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From this point on, we write a _L b for a , b relatively prime. We define the 
set of Ford circles: 

V ■— {C a ,b : a, b £ Z, a 1 b}. 

We define a set of circles to be normal if each circle is normal and no two 
circles have intersecting interiors (so, of course, V is normal). We order the set 
of normal sets of circles: 

A<B^ (J Ac U 5 ' 

AgA bgb 

Proposition 2.1. V is maximal with respect to the order <. 

Proof. We first show that no circle in V can be enlarged in the sense of having 
its radius increased but its tangent point to the a:-axis the same. Note that any 
two Ford circles have disjoint interiors (since \ad — bc\ > 1 for two Ford circles 
C a ,b and C C 'd)- Given a Ford circle C a< b, a _L b and so there exist c,d(d> 0) 
such that \ad — bc\ = 1 and thus there exists another Ford circle, C Cj( i, tangent 
to C a ,b- This implies that no circle in V can be enlarged and still maintain 
normality (i.e., its interior is still disjoint from the interiors of all the other 
circles in V). 

Next, we show that no normal circle can be added. Suppose there exist 
x,r so that for all Ford circles C a ,b, C(x,r)° D C a ^ = 0. Obviously x must be 
irrational and thus nx mod 1 is dense in [0,1], Consequently, the set {bx — a : 
a, b £ Z, b > 0} is dense in R and so, for any fixed r, there exist a, b such that 

a \/2r 

b < ~~b~~ 

By equation (1), this implies C(x, r)° (~l C° b ^ 0, a contradiction. □ 

We shall use the following geometric fact. 

Proposition 2.2. Given two tangent normal circles (tangent to the x-axis at, 
say, x and y), there is a unique third normal circle tangent to both and to the 
x-axis at a point between x and y. 

Proof. Suppose C(x,r)\\C(y, s) where x < y. For t £ {x,y), let Ri(t) := (t — 
x) 2 /Ar and R^t) = ( t—y) 2 /As. Then C(t, Ri(t))\\C(x, r) and C(t, R 2 (t))\\C(y, s) 
(and each is unique). As t increases, i?i(f) increases from 0 and i? 2 (i) decreases 
to 0 and so there exists a unique z £ (x,y) such that Ri(z) = R. 2 (z). Hence 
C(z,Ri(z)) is tangent to both C(x,r) and C(y,s). □ 

We call this new circle the child of the other two circles. It is easy to see 
that if C a ,b\\C c ,d then C a +c,b+d and C a - c ,b-d are tangent to both C a ,b and C c ,d- 
We say that C a ,b and C Ct d are parents of C a + c ,b+d and C a - c ,b-d 

It is possible to find the parents of a given Ford circle. For an ordered pair 
(a, b) of positive integers , consider the slow Euclidean algorithm: 

! [a, b — a] if a < b, 

[a — b, b] if a > b , 
stop if a = b. 
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This algorithm must terminate (since the sum of the two entries is positive and 
strictly decreasing) and, since the greatest common divisor is preserved at each 
step, this algorithm terminates with [gcd(a, 6), gcd(a, 6)]. For example: 

[14,5] A[9,5] A[4, 5] A[4,1]-A[3,1]^[2,1] A[l, 1]; (2) 

here we labeled each arrow according to which of the two entries is changed. 

This leads to a definition of maps L(a , b) := (a — b , b) and R(a, b) := (a, b — a) 
so that, for example, L 0 L 0 L 0 R 0 L 0 L(14, 5) = (1,1). The maps L and R 
are invertible and thus L _1 o L -1 o i? -1 o L -1 o i _1 o L _1 (l, 1) = (14, 5). 

In general, every relatively prime pair of positive integers (a, b) gives rise 
to a word w\W 2 ---w n in { L,R }* so that wf 1 o wf 1 o ... o tu“ 1 (l,l) = ( a,b ). 
It follows easily by induction that if (x,y) := wf 1 o wf 1 ° ° tu^ 1 (0,1) and 

(u, v) := wf 1 o W 2 1 o ... o rc“ 1 (l, 0), then C x , y and C UlV are the parents of C a ,b- 
We call this the parent algorithm. 

We now give a geometric construction of a normal set of circles that turns 
out to coincide with V. li A and B are normal sets of circles, we write A<B 
if B = A U {C} where C is the child of two circles in A (as constructed in 
Proposition 2.2). Starting with the normal set of circles Q 0 := {C n ^i : n £ Z}, 
let Q be the union of all normal sets of circles Q' that are maximal elements of 
a finite chain Q 0 < • • • < Q'. 

Lemma 2.3. Q = V. 

Proof. Note that C- a .b has parents C- x y and C- UtV if C at b has parents C x>y 
and C UjV . Hence the existence of the parent algorithm shows that any Ford 
circle C a ,b is either in (if b = 1) or C a ,b has two “parents” (if b ^ 1) Hence, 
by an induction argument, every C a ,b is in Q and thus, by the maximality of V, 
the theorem follows. □ 


We introduce a new way to parameterize normal circles. Given a, b £ 
(with a + b > 0), let 


(a, b) := C 


1 


a + b' 2(a + b) J 

It is easy to verify that every normal circle is represented in the form (a, b): 

1 — X X 


C{x,r) = 


2 r ’ 2 r 


and therefore every normal circle is represented in this new way. We refer to such 
representations as “barycentric” and will explain this terminology in Section 3. 

The set of Ford circles can be represented barycentrically in terms of solu¬ 
tions of a certain Diophantine equation. Let 

B := {(s, t) : s,t,u £ Z,gcd(s, t,u) = 1 ,s + t > 0, (s + t + u) 2 = s 2 + f 2 + u 2 }. 


Theorem 2.4. B = V = G■ 
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Proof. Given a Ford circle C a ^, let s := b 2 — ab, t := ab , and u := a 2 — ab. Note 
that (s,t) = C a ,b and it is easy to verify that (s,t) G B. Hence V < B and, by 
the maximality of V, and Lemma 2.3, the theorem follows. □ 

A surprising result follows (see [Ml). 

Corollary 2.5. If(a,b,c) is a relatively prime integer solution ofa 2 +b 2 + c 2 = 
(a + b + c) 2 then \a + b\ is a perfect square. 

Proof. If (a, b, c) satisfies the hypothesis, then (a, b) € B. By Theorem 2.4, 
(a, b) = C m>n for some m, n. The radii of these two circles are equal and 
therefore \a + b\ = n 2 . □ 

The results of this Section, with different proofs, have appeared in a paper 
[14] by the author and one of his students. 


3 Spheres, in general 


We identify R 3 with Cxi. For z £ C and r > 0, let S(z, r) be the sphere with 
center (z, r) and radius r. This sphere can be visualized as the sphere above and 
tangent to the complex plane at z with radius r. We say that a sphere is normal 
if it is of the form S(z,r ) for some z E C and r > 0. As in equation (1), it is 
easy to verify that S(z,r) and S(w,s) are tangent (we write S(z, r)||S'(w, s)) if 
and only if 

| z — w\ 2 = 47's. (3) 

Given three points in C, it turns out that there is a unique set of three 
mutually tangent normal spheres tangent to C at those points. 

Proposition 3.1. S(Pi,?r), <S(P 2 >r 2 ), and S(P 3 ,r 3 ) are mutually tangent if 
and only if 

|Pi-P 2 HPi~P 3 | I Pi - Pal • |Pa - P 3 \ IP1-P3MP2-P3I 
ri 2|P 2 -P 3 | ,r2 2|Pi — P 3 I ,r3 2|Pi — P 2 | 

Consequently, given three points Pi,P 2 ,P 3 , there exist unique real numbers 
r h r 2 ,r 3 > 0 such that S(Pi,r\), 5'(P 2 ,r 2 ), and S(P 3 ,rs) are mutually tangent. 

Proof. Suppose the spheres S(Pi,ri), S(P -2 , r 2 ), and S(P 3 , r 3 ) are mutually tan¬ 
gent. By (3), |Pi - P 2 | 2 = 4rir 2 , |Pi - P 3 | 2 = 4nr 3 , and |P 2 - P 3 | 2 = 4r 2 r 3 . 
Hence |Pi — P 2 | 2 /|Pi — P 3 | 2 = r 2 /r 3 and so 


\Pi 



P3 1 2 


— • 4r 2 r 3 = 4^. 

T3 


It follows that r 2 = 


|Pi - P2I ■ |P 2 - P 3 


and r 3 = 


2 |Pi-P 3 | 
Pi - P 3 | ■ |Pa ^ P 3 | 
2IP1-P2I 


. Similarly, r\ = 


I Pi ~ Pal • | Pi ~ P3 
2|P 2 - Pal 
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Conversely, suppose n = - 

\Pi - P 2 \ ■ | P 2 - P 3 


and r 3 = 


Pi-P 2 \- \Pi - P 3 1 
2\P 2 -P 3 \ 

Then 4riT2 = 


I Pi - Po\ ■ IA - P 3 


r 2 = 


IP - P 


2|Pi-P 3 | 

| 2 and so, by (3), 


21 _ p ^ | 

S , (Pi,r 1 )||S , (P 2 ,r 2 ). Similarly, 5(P 1; n)||5(P 3 ,r 3 ) and S(P 3 ,r 3 )||S'(P 2 ,r 2 ). □ 

It is well known that Mobius transformations take circles to circles (straight 
lines are considered circles too since, on the Riemann sphere, they are circles 
through oo; see m)- We shall use the notation 


(*) “ 


Every Mobius transformation C —> C extends to a unique sphere-preserving 
continuous map C x [0, oo) —> C x [0, oo) known as its Poincare extension. 
Its formula is available in ((’. for example. We shall use the fact that such 
an extension exists (even without needing its explicit formula) to express it in 
terms of normal spheres. Given a Mobius transformation m(z), extend m to 
normal spheres by defining 


: S(z 


S(m(z), |rn'( 2 :)|r). 


(4) 


Since 


| m(z) — m(w) | 2 = \m' (z)\\m'{w)\\z — w | 2 , 

it follows that if S(z, r)||S(w, s) then m(S(z,r))\\ih(S(w, s)). That is, rh pre¬ 
serves tangency of spheres. 

The Poincare extension rn of to, since it takes normal spheres to normal 
spheres, also acts as a function on the set of normal spheres: for some r' de¬ 
pending on z and r, 

m(S(z , r )) = S(m(z),r'). 

Given (w,t) £ C x R + , (w,t) is a point of intersection of two normal spheres, 
say S(z\,ri) and S(z 2 ,r 2 )- Let S(z 3 , r 3 ) be an arbitrarily chosen normal sphere 
tangent to both S(zi,ri) and S(z 2 ,r 2 )- Since {m(S(zi,ri)) : i = 1,2,3} and 
{rh(S(zi,n)) : i = 1,2,3} have the same three points of tangency to C, Propo¬ 
sition 3.1 implies m(S(zi,ri)) = m(S{zi,ri)) for * = 1,2, 3 and therefore 

rh(S(zi,n)) (~l rh{S{z 2 ,r 2 )) = {m(w , t)}. 

That is, the Poincare extension to of to is actually defined by (4). 

We now describe two ways to parameterize normal spheres. First, for a, fi £ 
C, let 

c ._ 1 


./3’2|/3| 2 , 

Clearly, every normal sphere can be uniquely represented thus: 

S{z,r) = 

We interpret Sqo to be the plane parallel, and one unit above, C. The Poincare 
extension is easily seen to have a simple formulation in terms of these spheres: 
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Proposition 3.2. If A := yj\ad — be\ then 


a b 
c 


Proof. Since 


• $(aa+bf3) / A,(ca+dfi) /A- 


a b\ : ^ ad — be 
c d) {Z) = (cz + d ) 2 ’ 


S, 


act+6y0)/A,(ca+d/3) / A 


= s 


aa + 6/3 |ad — 6c| 


ca + df 3 ’ 2|ca + d/3| 2 


= S' 


a b\ l a 
c d 


a 


1 


a b 
c d 


SI % 1 


2|/3p 


2|/3| 2 
a b 




□ 


A fact that will be used in the next Section is the following. 

Proposition 3.3. Given any three mutually tangent normal spheres, there are 
exactly two normal spheres so that each one, with the original three, form four 
mutually tangent spheres (with the understanding that planes parallel to the 
complex plane are spheres tangent at oo). 

Proof. Recall u := (—l + i\/3)/2 is a cube root of unity and that 0,1,1+w form 
the vertices of an equilateral triangle. The spheres So.i, Sip, Si+ W; i form three 
mutually tangent spheres each of radius 1/2. Clearly the plane z = 1 parallel to 
the complex plane (also denoted Si,o) is tangent to all the original three. The 
locus of tangent points of spheres that are tangent to two spheres chosen from 
So,i, Sip, Si+ 0,,1 form a straight line and so there is a unique sphere tangent to 
C at the intersection of three straight lines that is tangent to Sop, Sip, Si+ W p. 
Therefore the proposition holds for the three spheres. 

Now, suppose we have three mutually tangent spheres tangent at, say, Zi,Z 2 , Z 3 
respectively. There is a Mobius transformation taking these three points to 
0,1,1 + oj respectively and so, by Proposition 3.1, the Poincare extension takes 
the spheres to Sop, Sip, Si+ W p. Any sphere tangent to all three also get mapped 
to either S 1.0 or Sip_ w and therefore the inverse map applied to Si ,0 and Sip_^ 
gives all of the spheres tangent to those at z\, Z 2 , Z 3 ■ □ 

Recall the concept of “barycentric coordinates”. Given three non-colinear 
points in the plane (we henceforth define Po = 0, Pi = 1 ,P 2 = 1 + w so that 
the three points form an equilateral triangle of side length 1), it is possible to 
express every point in the plane uniquely as a real convex combination of the 
three: z = aPo+hPi+cP 2 , a+b+c =1. In this case, we say that z has barycentric 
coordinates (a, b, c). Removing the condition a + b + c = 1, we say that z has 
“projective barycentric coordinates” (a, b, c) if z = {aPo + bP\ + CP 2 )/{a + b-\- c). 



It is then possible to describe every normal sphere uniquely in terms of terms 
of three real numbers a, b , c where a + b + c > 0: 


(a, b, c) := S 


/ aPi + bP 2 + cP 3 
\ a + b + c 


1 

2(a T b -\~ 



We refer this parameterization as “barycentric” and note that every normal 
sphere has a barycentric representation: for real x, y, r with r > 0, 


S(x + iy , r) 


1 -x-y/Vi x-y/y/3 2y / y/3 ' 


2 r 


2 r 


2 r 


The barycentric parameterization allows a test for tangency. Consider the 
bilinear form 

Q(u,v) := (u- l)(v 1) - u v 
where u,v£l 4 and 1 := (1,1,1,1). Let 

U := {(a, b,c,d) el 4 : (a + b + c + d) 2 = a 2 + b 2 + c 2 + d 2 , a + b + c > 0}, 
and, for u := (a, b, c, d), let (u) be the sphere (a, b , c). 

Theorem 3.4. For u, v gM, 


Q{ u, u) = 0, 

Q(u, v) = 1 if and only i/(u)||(v). 

Proof. Since Q((a, b, c, d), (a, b , c, d)) := (a + b + c + d) 2 — ( a 2 + b 2 + c 2 + d 2 ), 
the first claim is obvious. 

Recall Pj (j = 1,2,3) form the vertices of an equilateral triangle with side 
length 1. Note that since |P ;| 2 + |Pj | 2 — 2SR(PjPj) = |Pj — Pj | 2 = 1 if i j, the 
norm squared of an R-linear combination of the Pi is 

|zPi + yP 3 + 2P3I 2 = (xPi + yP 2 + zP 3 ){xP 1 + yP 2 + zP 3 ) 

= £ 2 |Pi | 2 + y 2 \P 2 \ 2 + z 2 \P 3 \ 2 + 2xy$t(P 1 P 2 ) + 2xz^(P l P 3 ) + 2 yz$t(P 2 P 3 ) 

= (x + y + z)(x\P\\ 2 + y\P 2 \ 2 + 2 IP 3 I 2 ) - (xy + xz + yz). 

Hence, assuming a + b + c= l = s + t + u, 

|(aPi + bP 2 + CP 3 ) — (sPi + tP 2 + UP 3 ) | 2 

= |(a - 8)Pi + (b- t)P 2 + {c- u)P 3 \ 2 

= —[(a — s)(b — t) + (a — s)(c — u) + (b — t)(c — u)] 

= (a 2 + b 2 + c 2 + s 2 + t 2 + u 2 — 2 as — 2 bt — 2cu)/2. 

Let (a, b, c) and (s, t, u) be two spheres. Then the square of the distance 
between their tangent points is 

1 / a 2 + b 2 + c 2 s 2 + t 2 + u 2 2 (as + bt + cu) \ 

2 \ (a + b + c ) 2 (s + t + u) 2 (a + b + c)(s + t + u) J 
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and therefore they are tangent if and only if 

2 I j 2 I 2 2 I J.2 I 2 

, , . a + o + c , , , s +t +u 

(s + t + u) - - - 1 - (a + b+ c) -= 2(1 + as + bt + cu). (5) 

a b c s -\-1 -\- u 

Let d := —(ab + ac + be)/(a + b + c) and v := — (st + su + tu)/(s + t + u). 
Since (a + b + c) 2 = a 2 + b 2 + c 2 + 2 (ab + ac + be), it follows that a + b + c = 
( a 2 + b 2 + c 2 )/(a + b + c) — 2d and a similar result holds for s, t, u and v. The 
left hand side of equation (5) then becomes (s + t + u)(a + b + c + 2d) + (a + 
b + c)(s + t + u + 2v) = 2 (a + b + c + d)(s + t + u + v) — 2 dv and so the spheres 
are tangent if and only if 

{a + b + c + d)(s + t + u + v) = 1 + as + bt + cu + dv. ( 6 ) 

□ 

Remark. A very direct proof of this theorem is as follows. Given a := x + yui, 
/? := u+viv, 7 := X+Yu, and 6 := U+Vuj for real numbers x, y, u, v, X, Y, U, V, 

S a ,0 = (a,b, c) 

where a := u 2 + v 2 — uv + xv — xu — yv , b := xu — yu + yv, c:= yu — xv, and 
d := x 2 + y 2 — xy + xv — xu — yv. Similarly, S^ t g = ( A,B,C) for similarly defined 
A, B, C, and D. It is then an elementary but tedious calculation to see that 

M - /3t| 2 = Q{{a, b, c, d), ( A , B, C, D)). 


4 Ford spheres: tetrahedral case 



Figure 2: Ford spheres in tetrahedral arrangement. 

Just as Ford circles V were parameterized by Q, we shall parameterize a 
certain class of normal spheres V u by the Eisenstein rationals Q(w). Further, 
we shall show that the spheres of V u can be parametrized barycentrically by 
solutions of a certain Diophantine equation. A recursive geometric construction 
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of these spheres is also shown. In this way, we show an analogue of Theorem 
2.4 equating three parameterizations P u , B u , and Q u . 

After defining them appropriately, there are many possible strategies to 
showing they are all equal. Ours will be to show Q u C V w , V u C and, 
finally, B u C Q u . These spheres have been studied previously by Hellegouch 
and Rieger in [5].[ID]. [21] , and [25] , 

Recall oj := (—1 + iv3)/2 so that w 3 = 1 and w 2 = U. It is well known that 
the set of Eisenstein integers Z[w\ := {a + buj : a, b £ Z} is a Euclidean ring 
(see m)- norm is N(a + buj) = \a + bu | 2 = a 2 — ab + b 2 . As Z[w] is a UFD 
containing Z, not every prime in Z is prime in Z[oj]. It turns out though that 
p £ Z[ui] is prime if and only if either N(p) is the square of a prime congruent 
to 1 modulo 3 or equals a prime not congruent to 1 modulo 3 (see [11]). The 
Eisenstein rationals are the members of the field 


Q(w) := (r + soj : r, s £ Q} = {a/0 : a, 0 £ Z[w\,0 ^ 0}. 


Z[ui\ has exactly six units: ±1, ±ui, and ±ZU. We say that Eisenstein integers 
a, 0 are relatively prime (we write a T 0) if the only Eisenstein integers dividing 
both a and 0 are units. Since Z[uj\ is a unique factorization domain (UFD), 
a T 0 if and only if there exist two Eisenstein integers x, y such that ax + 0y 
is a unit. 

Recall that for a, 0 £ C, 



We define the set of Ford spheres: 


Vu ■■= {S a ,p : a,0 £ Z[ui\,a _L 0}. 


By Proposition 3.3, given three mutually tangent spheres, there is a fourth 
sphere so that these four are mutually tangent. We say that these four spheres 
are in a tetrahedral arrangement (since their contact graph is a tetrahedron). 
In fact, there are two ways to do this; we rephrase Proposition 3.3 thus: 

Lemma 4.1. There are exactly two tetrahedral arrangements of spheres con¬ 
taining three given mutually tangent spheres. 

We now define Q u recursively. Start with the set So := {So,i, ‘S'i+w,i}> 
perform the following process. Given S ni choose from it three mutually tangent 
spheres, and add the two as defined in Lemma 4.1 and thus form In 

this way, we have a “chain” of sets So C <Si C .... Let be the union of all 
elements in all such chains. Let the rank of a sphere S in denote the length 
of the shortest (finite) chain containing S. For example, rank(S'i i i) = 0 and 
rank(S'i j i_ w ) = 1. Hence a sphere of rank n has three “parents”, each of rank 
strictly less than n, so that all four are in a tetrahedral arrangement. Figure 2 
illustrates (some of) the spheres of rank at most 2: along with the three spheres 
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in So, we include the “child” Si t i- U of those three (the other child, Si t o, being 
a plane parallel to the complex plane, is not shown) as well as children of each 
of the triples formed by pairs chosen from <So and the sphere 5i ) i_ w . 

Lemma 4.2. C V u . 

Proof. Note that every element of of rank 0 is in V w . Suppose that every 
element of Q u of rank less than n is in V u and suppose S 6 Gu has rank n. Then 
its parents have smaller rank and are thus in V 0J . We will show that S £ V u . 

Let U := {1,1 + w,w, —1, —1 — w, —w} the set of all six units in Z[w]. Given 
three mutually tangent spheres S a> p, S 7 j, S XtV £ V , note that |a<5 — /3y| = 1. 
Letting w\ := x — a and w 2 := y — 0, the tangency of S x>y with the other two 
spheres implies 

aw 2 — 0W\ £ U, Sw\ — 7W2 = 0. 

Hence for some p, Wi = py and W 2 = pS. It follows that p(aS — /3y) is a unit 
and thus p is too. 

There are exactly two choices of a £ U such that cr — p £ U. Note that 
for any such a, S a + ait / 3 + a s together with the original three spheres forms a 
tetrahedral arrangement. Hence S = Sa+a-y^+a-s for some unit a and therefore 
S£V. □ 

Consider now the equation 

(a + b + c+ df = a 2 + b 2 + c 2 + d 2 . (7) 

We define l\ to be the set of relatively prime integer quadruples and we define 
a set of Ford spheres “barycentrically”: 

Bu := {( a,b,c ) : ( a,b,c,d ) £ (a + b+c+d,) 2 = a 2 +b 2 + c 2 + d 2 , a + b+c > 0}. 

Lemma 4.3. V u C B u . 

Pi'oof. Let S a jj £ V u . Then a 1 (3 where a = x + yw and 0 = u + vui for some 
x, y,u,v £ Z. Define 

a = u 2 + v 2 — uv + xv — xu — yv, 
b = xu — yu + yv, 
c = yu — xv, 

d = x 2 + y 2 — xy + xv — xu — yv. 


It is easy to verify 


Then 


(u + b + c + d) — ft T b~ c T d . 
a b + c + clo 1 1 


0 a + b + c 2|/3| 2 2(a + 6 + c)’ 

and so S a ,/3 = (a, b, c). 
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Since Z[u>\ is a principal ideal domain and a _L f3 in it, there exist 7 , 1 5 € Z[w] 
such that \a8 — /? 7 | = 1 and so S , a ,/ 3 || 5 ' 7 , 5 - 

As above, for some integer quadruple (A, B,C, D) satisfying (7), = 

(A,B,C). Hence (a,b,c)\\{A,B,C). By Theorem 3.4, 

(a + b + c + d){A + B + C + D) — (aA + bB + cC + dD ) = 1, 

and thus a, 6 , c, d are relatively prime. That is, (a, b , c) G B u and the lemma is 
shown. □ 


Theorem 4.4. B u ='P U1 = Q u . 

Proof. It is enough to show B :jJ C Q u . Consider the “generalized slow Euclidean 
algorithm” (‘gSEA’): 


(a, b, c, d) 


(—a, a + b, a + c,a + d) 
(a + 6 , —b, b + c,b + d) 
(a + c, b + c, — c, c + d) 
(a + d,b + d,c + d, —d) 


if a = min{a, 6 , c, d}, 
if b = min{a, b, c, d} < a, 
if c = min{a, 6 , c, d} < a, b, 
if d = min{a, b, c, d} < a, b , c. 


First, notice every solution of equation (7) is invariant under the transformation 
(a, b,c,d) 1 ^ (a + d,b + d, c + d, —d). Also, if d = min{a, b, c, d} then (a + d) + 
(6 + d) + (c+d) + (—d) = a + & + c + 2 d < a + b + c+d and, in general, a + b + c+d 
strictly decreases with every step of the gSEA. It cannot go below 0 however 
since a + b + c> a + b + c + d> 0 , 

a + 6 +c+ 2 d = a + b+c— 2(ab+ac+bc)/(a+b+c) = (a 2 + b 2 +c 2 )/(a+b+c) > 0. 

Therefore, if a + b + c+d > 0, then the gSEA eventually terminates with 
a, b,c,d > 0. If, say a, b > 0 at that stage, then 0 = ab + ac + ad + be + bd + 
cd > ab > 0 - a contradiction; therefore, at most one of a, b , c, d is non-zero. 
Since the ged is preserved by the gSEA, the gSEA eventually terminates at one 
of {g, 0,0,0), ( 0 , 5 ,0,0), (0,0, g,0), (0,0, 0,^) where g := gcd(a, b, c, d). Every 
u £ U must end in a basis vector ( 1 , 0 , 0 , 0 ), ( 0 , 1 , 0 , 0 ), ( 0 , 0 , 1 , 0 ) or ( 0 , 0 , 0 , 1 ) 
and reversing the gSEA in each of the other three basis vectors gives rise to 
three new parents since tangency is preserved by the gSEA. This is the “parent 
algorithm”. 

Given a relatively prime integer solution of (7), let its rank be the number 
of steps taken by the gSEA. For example, 

(12,12,3, — 8 ) ^ (4,4, -5, 8 ) ^ (-1,-1, 5, 3) ^ (1, -2,4, 2) 

^ (- 1 , 2 , 2 , 0 )^( 1 , 1 , 1 ,- 1 )^( 0 , 0 , 0 , 1 ). 

We code the steps: [4, 3,1, 2,1,4] and the rank of (12,12,3, — 8 ) is 6 . Reversing 
these (and noting that each step is idempotent), and applying to the standard 
basis vectors, we get: 

( 0 , 0 , 1 , 0 ) ^( 0 , 0 , 1 , 0 ) ^...^( 2 , 2 , 0 ,- 1 ) 

(0,1,0,0) ^(0,1,0,0) h4...^(5, 6 ,2, -4) 

(1,0,0,0) ^ (1,0,0,0) ^...^( 6 ,5,2, -4). 
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Note that the first step always leaves the vector fixed and so the parents 
(2, 2, 0,-1), (5, 6 , 2, —4), ( 6 , 5, 2, —4) of (12,12,3, — 8 ) have smaller rank. 

Note that every element of B u of rank 0 is in Q u . Suppose that every element 
of B u of rank less than n is in Q u and suppose S £ B u has rank n. Then its 
parents have smaller rank and are thus in Q u . It is easy to check that S € Q u . 
By induction on rank, the theorem is shown. □ 

Three corollaries immediately follow. 

Corollary 4.5. The integer solutions of 

(o + b + c + d) 2 = ci 2 T b~ -\- c 2 -f- d 2 


are parameterized by 


a = u 2 + v 2 — uv + xv — xu — yv , 
b = xu — yu + yv , 
c = yu — xv, 

d = x 2 + y 2 — xy + xv — xu — yv. 

The following Corollary solves a recent Monthly problem m 

Corollary 4.6. If a, b, c, d are relatively prime and satisfy (a + b + c + d) 2 = 
a 2 + b 2 + c 2 + d 2 , then |a + b + c| = m 2 + mn + n 2 for some integers to, n. 

Proof. Without loss of generality, assume a + b + c > 0. By hypothesis, (a, b, c) 
is a Ford sphere with radius l/2(a + 6 + c). By Theorem 4.4 , this sphere is also 
of the form S a: p with radius 1/2|/3| 2 for some Eisenstein integers a, ft. Hence 
\a + b + c\ = |^| 2 = to 2 + mn + n 2 for some integers to, n. □ 

Theorems 3.3 and 4.4 provide a useful way to construct Ford spheres recur¬ 
sively. 

Corollary 4.7. (Tetrahedral Rule). For a, b, c,d € U , if (a), (b), (c), (d) 

are mutually tangent, then so are (a), (b), (c), (a + b + c — d). 

Proof. Q(a, a+b + c-d) = Q(a, a) + Q(a, b) + Q(a, c)-Q(a, d) = 0 + 1 + 1-1 = 

1. □ 

The reason we call it the “tetrahedral rule” should be clear from Figure 3 
which represents several spheres as triples of integers obeying the tetrahedral 
rule. It is clear that this diagram indicates a tesselation of (part of) R 3 by 
tetrahedra such that every edge is shared by exactly six tetrahedra. The group 
of symmetries here is then a reflection group with four generators (Coxeter 
diagram a tetrahedron with each edge weight three (( ab ) 3 = e, etc.). This 
group is represented as a group of matrices Mk defined by 

(d//,- )ij = Sij P Sik SSikdjk- 
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010 


001 



421 


412 


Figure 3: Tetrahedral rule. 


This group seems related to the Eisenstein-Picard modular group [3]. 

An interesting side question is “when does the gSEA eventually repeat?”. 
The SEA of Section 2, starting with (1,#) eventually repeats if and only if the 
continued fraction for x eventually repeats (via the bijection (a, b) -o- b/a). 

It turns out that this is contained in the gSEA case. Note that for any x £ R, 
(1, x, x 2 , — x) is a solution of equation (7). We shall show that the gSEA applied 
to this vector eventually repeats if and only if x is a quadratic surd. To help, 
we say that two vectors A , B are equivalent (A = B) if a permutation of one is 
a scalar multiple of the other (e.g. (4,1, 3, 2) = (2,4, 6 , 8 )). 

Suppose first that x > 1. Then 

(1, x, x 2 , -x) (1 - x, 0, x 2 - x, x) (x - 1,1 - x, (x - l) 2 ,1) = (1, y , y 2 , -y) 

where y = x — 1 . 

On the other hand, if 0 < x < 1, then 


( 1 , x, x 2 , —x) i ^ ((1 — x, 0 , x 2 — x, x) 

^ ((1 - x) 2 , -x(l - x), x(l - x), x 2 ) = ( 1 , y, y 2 , -y) 


where y = x/(l — x). Hence the double-stepped gSEA is, modulo =, equivalent 
to iteration of 



This is closely related to continued fractions: if x has continued fraction expan¬ 
sion [ao, ai, 02 , 03 ,...], then it is easy to verify that 


fa 0 {x) = l/[ai,a 2 ,a 3 ,...] 
fao+ai (*^) [ 02 , O 3 , O 4 , ...] 

fa 0 +ai+a 2 {x) = l/[d3, O 4 , 05 , ...] 

etc. 

Hence the gSEA beginning with (1, x, x 2 , — x) eventually repeats if and only if 
x is a quadratic surd. 
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5 Ford spheres: octahedral case 



Figure 4: Octahedral Spherical Array 

In Section 4, we introduced an array V u of spheres parameterized by the 
Eisenstein rationals Q(w), a field whose ring of integers is a Euclidean domain. 
In this Section, we introduce an array V, of normal spheres parameterized by 
the Gaussian rationals Q(i), another field whose ring of integers is a Euclidean 
domain. In Section 4, we gave a recursive geometric procedure for constructing 
V u based on repeatedly, for any three mutually tangent spheres, adding another 
so that all four have contact graph a tetrahedron. In this Section, we shall give 
a recursive geometric procedure for constructing Vi that repeatedly, for any 
three mutually tangent spheres, adds three others, so that all six have contact 
graph an octahedron. Figure 4 shows three mutually tangent spheres together 
with their three children so as to form an octahedral array. In Section 4, we 
also gave a barycentric parameterization of the spheres in V u by expressing 
them in the form (a, &, c) for each relatively prime integer solution (a, b, c, d) of 
(a + b + c + d) 2 = a 2 + b 2 + c 2 + d 2 . By replacing this equation by 

(a + b + c + d) 2 = 2 (a 2 + b 2 + c 2 + d 2 ), ( 8 ) 

we get a new array of spheres, which under the Poincare extension M of a 
certain Mobius transformation, is the same as Vi. 

We shall first define Vi and, by using the (slow) Euclidean algorithm for 
Gaussian integers, show that it is maximal in the sense that no normal sphere 
can be included or enlarged. Next we show that octahedral arrays of spheres 
exist and that there are exactly two such arrays containing a given triple of 
mutually tangent spheres. This allows for a way of defining a set of spheres Qi 
geometrically in terms of octahedra just as we did earlier for tetrahedra. We 
then show that Vi C Qi (and so they are equal by the maximality of Vi). Next we 
define a set of spheres Bi, parameterized by solutions of (8), and show that V t C 
M{Bi ) where M is the Poincare extension of a certain Mobius transformation. 
We may thus conclude that Vi = Qi = M{Bi). 
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These spheres, parameterized as in the definition of Vi , have been studied 
previously by Pickover [231 . 

It is well known that the set of Gaussian integers Z[i\ := {a + ib : a,b G Z} 
is a Euclidean ring (see ED- Its norm is N(a + ib) = |a + ib\ 2 = a 2 + b 2 . As 
Z[i\ is a UFD containing Z, not every prime in Z is prime in Z[*]. The Gaussian 
nationals, are the members of the field 


Q(*) := {r + is : r, s G Q} = {a//3 : a,(3 G Z[i\,(3 ^ 0}. 


Z[i\ has four units: 1, —1 ,i, and —i. We say that Gaussian integers a, (3 are 
relatively prime (we write a T (3) if the only Gaussian integers dividing both a 
and /? are units. Since Z[i] is a unique factorization domain (UFD), a T f3 if 
and only if there exist two Gaussian integers x, y such that ax + (3y is a unit. 

Recall that for a, [3 G C, 



For this Section, we define the set of Ford spheres as 

Vi := {S a ,p : a, (3 G Z[i],a ± /?}. 


For any S a> p G Vi, a 1 f3 and thus there exist Gaussian integers 7 ,6 such 
that \a5 — /Jy | = 1. Hence S a ^\\Sry +pl ^+p5 for any p G Z^. For any two spheres 
G Vi, since aS — /3y is a Gaussian integer, \aS — /3y| > 1 and so 
S a ,p and S- f .s have disjoint interiors. Hence a given Ford sphere is tangent to 
infinitely many others but does not intersect the interior of any other. 

Lemma 5.1. For all z G C, 

\\\f{\f3z — a\ : a, f3 G Z[i]} = 0. 

Proof. We may define a “floor function” for C: 

[x + iy\ := [xj + i[y\ 

and a corresponding “fractional part” function 

{x + iy} := {x} + i{y}. 

Given positive integer N, divide the square R := [0,1) x [0,1) into N 2 disjoint 
congruent squares. For z ^ Q (i), the numbers in {( m+in)z :m,n= 1,..., iV+1} 
are all distinct and so, by the pigeonhole principle, there exists a, (3 G Z[i\ such 
that {az} and {f3z} are in the same small sub-square of R. There are thus 
7 , 6 G Z[i\ such that 


l(« - P)z + (7 - <5)1 = I {az} - {/3z}\ < 


Since N was arbitrary, the result follows for any z Q(i). The result obviously 


holds for z G Q(i). 


□ 
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Lemma 5.2. Vi is maximal. 


Proof. Note that any two Ford spheres S at p,S 7t $ have disjoint interiors (since 
\a5— /3y| > 1). Furthermore, since a _L /3, there exist 7 , <5 such that |a<5—/3y| = 1 
and so 5 aj) 8 || 5 7 i , 5 . Hence no Ford sphere can be enlarged. 

Suppose there exist z,r so that for all a,b, S(z,r)° fl S a $ = 0. Obviously 
2 ^ Q(*) and thus, by Lemma 5.1, the set {/3z — a : a € Z[w],/3 £ Z[w] + } is 
dense in C and so, for any fixed r, there exist a, /3 such that 

a y/2r 

1 < W 

This implies 5(z,r)°n5° / 3 ^0,a contradiction. □ 

We now define a set of spheres Ri recursively. Start with the set So := 
{ 50 . 1 , 51 , 0 }, perform the following process. Given S n . create 5„+i by choosing 
two tangent spheres, say 5 a> / 3 , 5 7l i in S n and add the sphere S a +py,p+p6 for 
one of the units p £ {1,-1 In this way, we have a “chain” of sets 

5 0 d S 1 (Z .... Let 'R., be the union of all sets in all such chains. 

Lemma 5.3. Vi = 7 Zi. 

Proof. By Lemma 5.2, It is enough to show that V, C 72.,. We shall do this by 
utilizing the “slow Euclidean algorithm” for the Gaussian integers. 

Given Gaussian integers a and /3, one can choose a unit p so that the angle 
between a and p/3, as vectors, is in [— 7 t/ 4 , 7 t/ 4 ], Therefore, if 0 < |/3| < |a|, 
then there is a unit p such that \a — p/3| < |a|. We now define SEAi, the slow 
Euclidean algorithm for Z [?’]: 

! [a — p/3,/3], if 0 < |/3| < |a| and p £ U with |a — p/3| < |a| 

[a,/3 — pa], if 0 < |a| < |/3| and p £U with |/3 — pa| < |/3| 
stop, if a/3 = 0. 

Note that |a | 2 + |/3 | 2 is decreasing, positive, integer valued so the algorithm must 
stop for any two Gaussian integers 01 , (3. Since the “greatest” common divisors 
are preserved at each step, the algorithm must end at either ( 0 , p) or (p, 0 ) for 
some p. Further, if a T /3, then p must be a unit. We say that the “rank” of a 
pair [a, /3] is the number of steps taken by SEAi to end. 

As with the parent algorithm for Ford circles, run SEAi until the penulti¬ 
mate state [pi,P 2 ] (where p\,pn € U) and then reverse the steps on [pi, 0 ] and 
[0, P 2 ]- The resulting pairs, call them [ai,/3i] and [a 2 ,^ 2 ] are “parents” of [a,/3] 
in the sense that each has lower rank than [a, /3] and the two must sum to [a, /3]. 

Every sphere of rank 0 is in Ri. Suppose that every sphere of rank less than 
n is in IZi and suppose 5 aj p £ Vi. Then its parents have lower rank and so are 
in 7 Zi. By the definition of IZi, S a ^p £ R,. By (strong) induction, the lemma is 
shown. □ 


z — 
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Given three non-colinear points Pi, P 2 , P 3 , we consider them on the bound¬ 
ary of a hyperbolic disk. There is then a unique geodesic triangle with vertices 
Pi, P 2 , P 3 and, a unique circle inscribed in that triangle. This gives rise to three 
unique points Qi, Q 2 , Q 3 such that the six points Pl, P 2 , P 3 , Qi, Q 2 , Q 3 form an 
octahedron (that we call a Mobius octahedron generated by Pi, P 2 , P 3 since a 
Mobius transformation changes it to another such octahedron) - see Figure 4. 



Figure 5: Mobius Octahedron 


The contact graph of a finite collection V of normal spheres is the graph with 
vertex set V where two vertices share an edge if the two spheres are tangent. 
For example, three mutually tangent spheres have tangency graph a triangle. 

Theorem 5.4. The six points of a Mobius octahedron are the tangent points of 
six spheres with contact graph an octahedron. 


Proof. Recall the cross ratio: 

[z,q,r, s] 


(z - q)(r - s ) 
[z - s)(r - q)' 


It is well-known that it is invariant under Mobius transformations: 


[m(z),m(q),m(r),m(s)\ = [z,q,r,s] 


(see Needham US!)- Let A, B , C, D , E, F be vertices of a Mobius octahedron 
where A,B,C form a triangle and D,E,F form a triangle. Let ( AB ) denote 
|A — B |, etc. It is easy to verify that 


\[F,C,B,D][A,B,D,E][A,C,F,E]\ 


(. AB)(AC){DE)(EF) 
(AE) 2 (BG)(DF) 


and so the right side is invariant under Mobius transformations. For the special 
case where A = 0, B = (i + l)/2, C = 1, D = 1 + i, E = 00 , and F = i, the 
right side is 1 and since every Mobius octahedron is the image of the special one 
under a Mobius transformation, 


(AE) 2 


(AB)(AC) {ED) (EF) 
(BC)(BF) 


(9) 
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By Proposition 3.1, there are three mutually tangent normal spheres tangent 
to the plane at A,B,C and another three tangent at D,E,F. By (9) and 
Proposition 4.1, |A — E \ 2 = Ataxe and so the spheres at A and E are tangent. 
Similar arguments for pairs AF, BD 1 BF , CD 1 and CE show that the six spheres 


have an octahedral tangency graph. 


□ 


“Extended” cross-ratios like those in the proof of Theorem 5.4 make an 
appearance in the paper 1111- 

Lemma 5.5. Given three mutually tangent spheres, there are exactly two ways 
to choose another set of three mutually tangent spheres so that those six have 
contact graph an octahedron. Further, every collection of six normal spheres 
with contact graph an octahedron are tangent to the plane at six points in some 
Mobius octahedron. 

Proof. Given three mutually tangent spheres, let Z\,Z 2 , Zs denote where they 
are tangent to C. The Mobius transformation m{z) := [z, Z\, Z 2 , Z 3 ] defined 
in terms of the cross ratio takes Zi, Z2, 23 to 0 , 1 , 00 respectively (see Figure 
6). Its Poincare extension takes the three spheres to S'0,1, Si,i, Si,o- It’s easy 
to verify that these three spheres together with the spheres S p>p +i, S Pt 1, <Si+p,i 
(where p = i or — 1 ) forms an octahedral arrangement and so the reverse Mobius 
transformation m~ l extends either one to an octahedral arrangement. There 
are no more than these two since the points 0 , 1 ,00 form a triangle in a Mobius 
octahedra if and only if that octahedron has vertices {0, 1, 00, i, 1 + i, (1 + i)/ 2} 
or {0,1,00, — i, 1 — i, (1 — f)/2}. 

Given six normal spheres with contact graph an octahedron, choose three 
mutually tangent spheres that are tangent to C at, say, z\,Z2,z^ . The image 
of all six spheres under the Poincare extension of [-, z\, Z2, Z3] is then, by the 
argument above, vertices of a Mobius octahedron (see Figure 6) and the result 
follows. □ 





Figure 6 : Image of Fig. 4 under [-, Pi, P 2 , P 3 ] 


We may thus define a new type of Ford sphere geometrically. We now define 
Qi recursively. Start with the set Sq := {<So,i, <S'i,o } 5 perform the following 

process. Given S n , choose three spheres, and add the three as defined in Lemma 
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5.1 and thus form In this way, we have a “chain” of sets So C C .... 

Let Qi be the union of all sets in all such chains. 

Lemma 5.6. 7 Zi C ft. 

Proof. We define, for U := (a, ft) £ Z[i] 2 , S[[7] := S a p. Note that for any unit 
p,S\pU] = S[U]. 

It is enough to show that if S[f7],S[V] are two tangent spheres in Qi and p 
is a unit, then S[U + pV] £ Qi. Suppose S[U], 6’[V] are two tangent spheres in 
Qi. The two spheres are then part of an octahedral arrangement and so part of 
a triangle which, without loss of generality, is of the form S[U], S[U + V], S[V]. 
It is easy to verify that for p = ±i, X p := {£[[/], S'[[/ + V], S[U], S[U + p(U + 
V)], S[U + pV],S[U + V + pV}} forms an octahedral arrangement of spheres and 
thus, by Lemma 5.5, these are the only octahedral extensions of S[f7],S[t/ + 
V\, S[U]. In this way, we see that S[U + pV) £ Qi for p = l,i, —i. The triangle 
of spheres S[U], S[U + iV\, S[U] is part of two octahedra, one of which contains 
S[U — V] and so the lemma is shown. □ 

We now reconsider equation ( 8 ): 

(a + b + c + d ) 2 = 2 (a 2 + b 2 + c 2 + d 2 ). 

A solution of ( 8 ) is known as a Descartes quadruple and it was known (by 
Descartes) that that the curvatures of four mutually tangent circles satisfy ( 8 ) 
(see [ 6 ]). We shall parameterize the sets set of spheres Vi by relatively prime 
integer solutions of ( 8 ). 

As before, we define Z” to be the set of relatively prime integer n-tuples. 
Let V denote the set of relatively prime solutions of ( 8 ): 

V := {(a, 6 , c, d) £ Z^ : (a + b + c + d) 2 = 2 (a 2 + b 2 + c 2 + d 2 )}. 

A related set, the set of “Descartes triples”, can be defined: 

S := {(a, b, c) e Zi : Vab + ac +be £ Z}. 

We call them Descartes triples since it is easy to see that if (a, 6 , c) £ S , then 
(a, b,c,a + b + c± 2 \/ab + ac + bc) £ V and, conversely, if (a, 6 , c,d) £ V then 
(a, b , c) £ S. 

Let 

M(z) := iz/(( 1 — z)oj) and M~~ x (z) := coz/(coz + i). 

Then M is the Mobius transformation taking the points 0,1,1 +uj (the vertices 
of an equilateral triangle of side length 1 ) to the points 0 , oo,i respectively. 
Using the notation of Section 3, 

Lemma 5.7. Given a, ft £ C, let a := \ft\ 2 +£s(aft),b := \a\ 2 +$s(aft),c := 
—$s(aft), and m = $t(aft). Then 

M~ 1 (S oli p) = (a + m/V 3, b + m/V 3, c + m/Vi) 

and \m\ = Vob + ac + be. 
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Proof. By Section 3, 

M-\S a , p ) = S ua , ua+if) = S f—-- 1 —^) • 

\tua+ ip 2\uja + ip\ z J 

For this sphere to agree with ( A , B,C) := S((B + C(1+uj))/(A+B + C), 1/2 (A+ 
B + C)), we must have 

A + B + C = \u>a + z/3| 2 = |a| 2 + |/3| 2 + v^»(a/3) + 3(a/3) 

and 

B + C( 1 + w) = toa/oja + i/3) = \a\ 2 — iuja/3 

= H 2 + \y/®t{7xp) + ^(a(3) + i Qsft(a/?) - i^33(a/3)) . 

Solving for A, B , C, we find 

C = $l(a(3)/V3 — 3 (a/3), 

13 = |a| 2 + $t(a/3)/y/3 + 5(a/3), 

A = |/3| 2 + 5ft(a/3)/V3 + 3(a/3), 

and the first part of the lemma follows. 

With a, 6, c defined above, it is easy to verify that 

ab + ac + bc— |a| 2 |/3| 2 — 3(a/3) 2 = 5i(a/3) 2 = to 2 . 

□ 

We now define a type of Ford sphere barycentrically: let 
Bi := {(a + m/V 3, b + m /V5, c + m/V3) : (a, b, c) £ S,m = ±\/ ab + ac + 6c}. 
Lemma 5.8. Pi C M(Bi). 

Proof. Define 

L(a, 6, c) := (a + m/Vs, b + m/Vs, c + m/Vs, — 2m/y[3) 

where to = \/ab + ac + be. It is easy to check that V (a, 6, c) satisfies equation 
(7). Further, with Q the bilinear form defined in Section 3, 

Q(V(a,b,c),V(A,B,C)) ^ 

= aB + aC + bA + bC + cA + cB — 2 \J ab + ac + bcV AB + AC + BC. 

Let S a ^i £ Vi. Then a _L /? £ Z[i\ and so |a| 2 , |/3| 2 ,5ft(a/3), 3(a/3) £ Z. By 
Lemma 5.7, 

S a ,p = M((a + m/V 3,6 + m/V 3, c + m/VH)) 

for integers a, 6 , c. Since a _L f3 there exist Gaussian integers 7 , <5 so that \aS — 
/3y| = 1 and also 7 _L S. Then S^ t s = M((A + n/V3,B + n/V3,C + n/V 3)) for 
integers A, 13, C, n. By (10), aB + aC + bA + bC + cA + cB — 2 mn = 1 and thus 
(a, b,c) £ S and thus S a> p G 

□ 
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Together, the lemmas of this Section imply 
Theorem 5.9. Gi = Vi = M(Bi). 

Corollary 5.10. For a Descartes quadruple ( a,b,c,d ), 

a + b + c = | 7 | 2 + |<5| 2 


for some 7 ,SG Z[w], and 

a + b = m 2 +n 2 

for some m, n £ Z. 

Proof. Using the proof above, 

a + b + c = |a | 2 + \/3\ 2 + 3(5/3) 

= (ft 2 + CI1&2 T ^2) “b ( a 2 — ®2^1 + b 2 ) = |fll — ^2 ^| 2 + |a2 + b\U}\ 2 . 

Given Descartes triple (a,b,c), let S(z,r) = M({a + m/y/3, b + m/y/3, c + 
m/y/ 3)) where M is as above and m 2 = ab + ac + be. Then 2 is a Gaussian 
rational and there exist relatively prime Gaussian integers a,j3 with z = a/(3. 
Then a = |/3| 2 + 3(a/3), 6 = |a | 2 + 3(a/?), and so 

a + b = |a | 2 + |/3 | 2 + 23(5/3) = (ai + 6 2 ) 2 + (a 2 - 61 ) 2 . 


□ 

The close relation between continued fractions and Ford circles has long 
been noted (see m for a thorough study of the relation between continued 
fractions and the dual of Ford circles). Higher dimensional analogues have also 
been studied by many authors (generally by extending the compact interval to 
a triangular region). Another option, not previously studied, is to extend the 
theory of continued fractions to the “Sierpinski gasket”, a natural generalization 
of the closed interval. In attempting this, one may try to generalize Ford circles 
to spheres on the Sierpinski gasket. 

A compact interval can be represented 


[Pi,P*\ 




2~ j Pf U) :f€{l,2y 


This generalizes easily to any dimension; in dimension 3, we have the “Sierpinski 
gasket” 


SG := 


:/€{!, 2 -3} N 


where Pi, P2 , P3 are complex numbers forming the vertices of an equilateral 
triangle. 
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Figure 7: Sierpinski Gasket SG and Apollonian Circle Packing CP 


The Sierpinski gasket is obviously homeomorphic to part of an Apollonian 
circle packing. The two figures in Figure 7 are called fractals because their 
Hausdorff dimension is fractional. Roughly speaking, since SG is made up of 
three copies of itself, with each copy having length and width half as big, the 
dimension of SG is the solution of 2 d = 3 (namely d = In3/In2). However, 
Hausdorff dimension is not a topological invariant and, indeed, the dimension 
of CP has been found to be 1.305688± 10~ 6 by McMullen [21] but is not known 
with complete precision (see Graham et al. [16]). 

A local cut point a; in SG or in CP is an element for which there exists a 
connected neighborhood U of x for which U — {ir} is disconnected. In SG, such 
points are of the form ]Tb 2where the corresponding function f(j) is 
constant for all sufficiently large j or, equivalently, the point can be represented 
by more than one /. 

The author’s original motivation for this study was the extension of contin¬ 
ued fractions on an interval (and their interpretation in terms of Ford circles) 
to a development of continued fractions on SG where, it was hoped, that a type 
of “Ford sphere” could be attached to each local cut point. It is indeed possible 
for CP. 

Given three mutually tangent circles C\ , C2 , C3 in C with respective tangency 
points {vjij} := Ci fl Cj, there is a unique set of three spheres S12, S 13, S23 that 
are tangent to C at W12, W13, W23 respectively. By Lemma 2 of [20], if Ci is the 
curvature Ci for i = 1,2,3, then the curvature of Sij equals Cj + Cj. 

By Theorem 5.4, it is clear that it is possible to assign a sphere to each 
local cut point of the Seirpinski gasket and, in general, to each tangency point 
of a “weak” circle packing such as an Apollonian super-packing (see 0)- The 
(0,0,1,1) super-packing is pictured below (from Figure 4 of [7]). 
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Figure 8: Part of the (0,0,1,1) super-packing 


6 Ford spheres: general case 

A complex quadratic number field is a field of the form F{D) := Q(y/—D). 
The set of algebraic integers of F(D) is 0(D ) := Op(D) = Z[cr] where a = 
(1 + y/—D)/2 if D = 3 (mod 4) or er = V ~D if D = 1,2 (mod 4) (see, for 
example, SB, S3])- It follows that 

l2 ( m 2 + mn + ^4—n 2 if D = 3 (mod 4) 

TO + TOT = < „ „ 

1 to 2 + Dn 2 if D = 1, 2 (mod 4) 


is, in all cases, an integer. 

From a celebrated theorem of Heegner, Stark, and Baker, for D > 0, O(D) 
is a principal ideal domain (equivalently, a unique factorization domain) if and 
only if D is one of the “Heegner numbers” {1, 2,3, 7,11,19,43,67,163}. Further, 
it is known that O(D) is a Euclidean domain (equivalent, in this case to norm- 
Euclidean) if and only if D is one of the first five Heegner numbers; i.e., D € 
{1,2,3,7,11} (see pT3] ). Gauss first posed the question, open to this day, of 
whether there are infinitely many real quadratic number fields (where D < 0) 
with class number 1. See, for example, m for a nice account of Euclidean 
domains and m for a more general account of quadratic number fields. 

For D e {1,2,3,7,11,19,43,67,163}, let er be such that 0{D) = Z[cr] and 
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define Ford spheres in this case to be 

V a := {S Q)J g : a,/3 £ Z[a],a _L /?} 

which thus forms an array of spheres with non-overlapping interiors and such 
that every sphere is tangent to many others. 

As for Ford circles and the Ford spheres of Sections 2 and 5, we shall show 
that V a is maximal. 

Lemma 6.1. For all z £ C, 

inf {| f3z — a\ : a, /? £ Z[ct]} = 0. 

Proof. We may then define a “floor function” for C: 

\x + <jy\ := L^J + 

and a corresponding “fractional part” function 

{x + ay} := {a:} + a{y}. 

Given positive integer N, divide the parallelogram R := {x + ya : x, y £ [0,1)} 
into N 2 disjoint congruent parallelograms similar to R. 

Clearly, the lemma holds for z £ Q(<r). Fix z £ C — Q(<r). For distinct 
integers j and k. {jz} ^ {kz} since, otherwise, (j — k)z € Z [a]. Hence, for j = 
l,...,n 2 + 1 , the numbers {jz} are distinct and so, by the pigeonhole principle, 
one of the little parallelograms contains distinct {jz} and {kz}. Hence, for some 
aj,a k £ Z[cr], 

|{jz — aj) — (kz — a k )\ = | {jz} — {kz }| < diam(R)/n 
and the lemma follows. □ 

Lemma 6.2. V a is maximal. 

Proof. Note that any two Ford spheres S a> p,S lt s have disjoint interiors (since 
\aS— /3y| > 1). Furthermore, since a _L f}, there exist 7 ,5 such that |a<5—/3y| = 1 
and so <S aj ^||5 7 >( s. Hence no Ford sphere can be enlarged. 

Suppose there exist z,r so that for all a, b, S(z,r)° D S ai p = 0. Obviously 
2 ^ Q(cr) and thus, by the previous lemma, the set {/3z—a : a £ Z [a\,/3 £ Z[cr] + } 
is dense in C and so, for any fixed r, there exist a, /? such that 

\/2r 

w 

This implies S(z,r)° fl 5°^ ^ 0, a contradiction. □ 
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A barycentric representation of these spheres are based on solutions of the 
equation: 

j ab + ac + bc+ (a + b + c)m. = ^=2 m 2 if D = 3 (mod 4) 

1 ab + ac + be = Dm 2 otherwise. 


( 11 ) 


Let B a := {(a + mf, b + m£, c + m£) : (a, b, c) £ Zj} is a solution of (11)} where 

(y/Z - y/D)/y/l2 if D = 3 (mod 4) 
y/~D /y/Z otherwise. 


£ = 


Let fi(z) := ^ ( 2 :) be the Mobius transformation taking the points 0,1, 00 

to 0,1 + w, 1 respectively. 

Theorem 6.3. n{P a ) = B a . 

To prove this, we first prove several lemmas. 

Lemma 6.4. //a/3 = s + it then 

h(Sa,i3) = (a + m/ V3, b + m/y/Z, c + m/y/Z) 
where a = |/3 | 2 — s, b = \a\ 2 — s, c = s, and ab + ac + be = m 2 = t 2 . 

Proof. Let a = |/3 | 2 — s, b = |a | 2 — s, c = s, and m = t. Then 

ab + ac + bc= |a/3 | 2 — s(|a | 2 + |/3 | 2 — 2s) = |a/3 | 2 — s 2 = t 2 = m 2 . 

Note that 

wa/3 — —s/2 T y/zt/2 H- i{y/Zs/2 -t- 1/ 2) 

and thus 

(a + m/ y/Z) + (b + m/ y/Z) + (c + m/ Vz) = |a | 2 + |/3 | 2 — s + ty/ 3 
= |a | 2 + |/3 | 2 + 25ft(a;a/3) = |wa + /3| 2 , 

and 

(b + m /Vz) + (c + m /V / 3)(1 + w) = |a | 2 - s/2 + v / 3t/2 + i(VZs/2 + i/2) 


= |a | 2 + wa/3 = wa(w« + /3). 
Hence, by Proposition 3.2 

/i(*S/. ) ,g) — S u ;a,wa+/5 — S 


W« 


= S' 


= s 


wa + /3’ 2|wa + /3| s 

wa(w« + /3) 1 \ 

|wa + /3 | 2 ’ 2 |wa + /3 | 2 J 

b + m/y/Z + (c + m/y/ 3)(1 + to) 


1 


y a + 6 + c + m-\/3 2(a + b + c + my/Z) 

= (a + m /a/3, b + m/y/Z , c + m /Vz). 


□ 
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Lemma 6.5. If fi(S a ,p) = (a + m/y/3, b + m/y/3, c + m/y/S) and /x(5 7 j( 5 ) = 

(a' + m! /V3, b' + m' /V3,c' + m' / \/3) t/ien 

|ad — /3y | 2 = ah' + ad + ba' + bd + co! + cb' — 2mm'. 

Proof. By Lemma 6.4, the hypothesis implies a = |/3 | 2 — s, b = \a\ 2 — s, c = s, 
m = t where a/3 = s + it and a' = |<5 | 2 — s', b' = |y | 2 — s', d = s', m' = t' where 
7<5 = s' + it'. Then 

a(b' + d)+b(a' + d) + c(a' + b') — 2mm' 

= m 2 - s)|7| 2 + (H 2 - s)|5| 2 + s(| 7 | 2 + |5| 2 - 2s') - 2W 
= \a5\ 2 + |/? 7| 2 - 2 (ss' + tt') 

= \a6\ 2 + |/? 7| 2 — 2$t(a/3'y6) = |cr<5 — /? 7 | 2 . 


□ 

Proof of Theorem 6.3. Suppose D = 3 (mod 4), a = (1 + iy/~D)/2. Given 
a, (3 £ C, let s + it = a/3, A := |/3 | 2 — s + t/y/D, B := \a\ 2 — s + t/y/D, 
C := s + t/y/D, M := —2t/y/D, and £ := (-\/3— y/D)/y/ 12. By Lemma 6.4, 

h(.d a ,p) = {A + Mf, B + All ;, C + Alf). 

By Lemma 6.5, if /z(5 a ,/ 3 ) = {A + M£, 33 + M£, C + M£) and y(S lt s) = (A' + 
M'£, B' + M'£, C' + M'f) then 

|aS - /3y | 2 = AB' + AC' + BA' + HC" + CA' + CB' 

n _ q ( 12 ) 

+ M (A' + B' + C') + M'(A + 33 + C)- —MM'. 

Given two distinct spheres in B a , say {A + M £, 33 + M£, C + M£) and (A' + 
M'f,B' + M'£,C" + M'£), their images under /i ” 1 are of the form S XtV ,S UtV 
respectively for some complex x,y,u,v. Equation (12) implies \xv — yu\ 2 is a 
(positive) integer and therefore the collection of spheres B a is normal (i.e., no 
two spheres have intersecting interiors). 

If a,/3 £ Z[er], then |a | 2 and |/3 | 2 are integers. Also, a/3 = s + it where 
2s,2t/y/D £ Z and 2s = 2 t/y/D (mod 2). Then A,B,C, and M are integers. 
If a -L /3 then there exist 7 , <5 such that |a<5 — /3y| = 1 and so, by (12), 

AB’ + AC' + BA' + BC' + CA' + CB' 

+ M{A! + B' + C') + M'{A + B + C) - MM' = 1 

for some integers A',B',C',M' and therefore A,B,C,M are relatively prime. 
Hence /z(S Q)l g) £ B a . By the maximality of V a , the theorem holds. 

Suppose D = 1,2 (mod 4), a = iy/D. Given a,/3 £ C, let A := |/3 | 2 — s, 
33 := |a | 2 — s, C := s, M := t/y/D, and f := y/D/y/3. By Lemma 6.4, 

h(.d a ,p) = {A + Mf, 33 + Ml ■, C + Al^)■ 
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By Lemma 6.5, if y(S a >( g) = (A + M£, B + M£, C + M£) and /j,(Sy t g) = {A' + 
M'£, B' + M'£, C' + M'£) then 

| aS - / 37| 2 = AB' + AC" + BA' + BC' + C'A' + CB' - 2D MM'. (13) 

Given two distinct spheres in B a , say (A + M£, 13 + Mf, C + M£) and (A' + 
M't;,B' + M'£,C' + M'£), their images under /x _1 are of the form S X}V ,S U , V 
respectively for some complex x,y,u,v. Equation (13) implies \xv — yu\ 2 is a 
(positive) integer and therefore the collection of spheres B a is normal (i.e., no 
two spheres have intersecting interiors). 

If a,/3 £ Z[er], then \a\ 2 and |/3 | 2 are integers. Also, a/3 = s + it where 
s,isZ. Then A, B 1 C, and A/ are integers. If a _L /3 then there exist 7 , S such 
that |aS — (3j\ = 1 and so, by (13), 

AB' + AC + BA' + BC' + CA’ + CB’ - 2D MM' = 1 


for some integers A! , B',C , M' and therefore A,B,C,M are relatively prime. 
Hence 3 ) G B a . By the maximally of V a , the theorem holds. □ 


As in the proof of Theorem 6.3, since a 1 /3 implies A, B , C, M are relatively 
prime, we have a test for relative primality in Z[er]. 

Corollary 6.6. a T p in Z[cr] if and only if 

gcd(|a| 2 , |/3| 2 , s,t/y/D) = 1 if D = 1,2 (mod 4) 

gcd(|a| 2 , |/3| 2 , s — tVD, 2t/V~D) = 1 if D = 3 (mod 4) 
where s + it = a/3. 


We now find an algorithm for finding the relatively prime integer solutions 
of (10) (and thus the elements of B a and V a ). For any quadratic polynomial 
f(x), the “secant addition” 


x © y 


xf(y) ~ yf{x) 

f(y) - /O) 


is associative and, in fact, if / has roots u and v, then ((M—{u, n})U{oo}, ®, 00 ) is 
an abelian group (see HZ] ). The ring of algebraic integers Z[cr\ has characteristic 
polynomial 


f(x) 


x 2 -x + ^- if D = 3 (mod 4) 
x 2 + D otherwise. 


The corresponding secant addition is defined by 


x © y := 


xy — (D + l)/4 

x + y — 1 
xy — D 

x + y 


if D = 3 (mod 4) 
otherwise. 


It follows that £©y© 2 i = ooif and only if (— x, —y, —z ) is a solution of (11). 
Hence, to find (all) ( a,b,c ) £ Z 3 ± that solves (11): 
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• Choose x, y £ Q, 

• Calculate z := (x © y) _1 , 

• Let m be the least positive integer so that m{xy + xz + yz) € Z, 

• Let a = —mx , b = —my , c = —mz. 

Just as there is a recursive geometric construction of B c , based on a tetra¬ 
hedron when D = 3 and based on an octahedron when D = 1, we believe 
that for all five cases D = 1,2, 3, 7,11 corresponding to Z[ct] being a Euclidean 
domain, there is a recursive construction based on a polyhedron. Figure 9 
shows the conjectured polyhedra with respective discriminants 3,4,7,8,11 (i.e., 
D = 3,1, 7, 2,11 respectively). Computer experimentation supports the conjec¬ 
tured relation, as do some results in the literature (see [32]). 



Figure 9: Polyhedra corresponding to the five imaginary quadratic Euclidean 
domains. 


7 Conclusion 

In this section we attempt to connect the material of this paper to past, present, 
and future research. 

The main direction of this paper was to understand and generalize Ford 
circles by, first, developing three different parameterizations of them {V,Q,B)\ 
the third being apparently new. There are therefore three different directions 
Ford circles to be generalized. We first consider the geometric, followed by the 
algebraic, and then “barycentric”. 

The initial motivation for this paper was to extend continued fractions in 
a new way. Since continued fractions are closely related to the geometry of 
Ford circles (a topic not developed in this paper but which appears, implicitly 
at least, in [3(5]), any geometric generalization of Ford circles could lead to a 
generalization of continued fractions. A natural extension of the interval is to 
the Sierpinski gasket (or any of its homeomorphic images) see Figure 7, and 
a generalization of continued fractions to that space might be new. Assigning 
spheres to the local cut points of CP seemed to be promising and indeed this 
is the type of sphere developed in Section 5. 

Problem 7.1. Develop a theory of continued fractions on CP. 
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Three-dimensional analogues of Ford circles have been studied previously by 
Hellegouarch mm, Pickover [23], and Rieger [24[ [25] but a unifying treatment 
of them has not been done. This paper, although not explicitly a survey of the 
topic, was meant to cover the topic more comprehensively. In particular, we 
cover the cases corresponding to Q(i), Q(w) and, to a lesser extent, Q(er) for 
which Z[cr] is a UFD. 

Problem 7.2. Develop the connection between continued fractions and the the 
family V a of spheres. 

Diophantine approximation is a major theme of continued fractions and can 
be understood geometrically in terms of Ford circles; see [3], p. 30. Schmidt 
[26l 27, [23 [29] developed the theory of Diophantine approximation for four 
of the imaginary quadratic number fields for which the ring of integers is a 
Euclidean domain. For a rational number q := a/b ( a,b in lowest terms), let 
C(q ) := C a ,b- The “parent algorithms” of Sections 2, 4, and 5 leads to an 
algorithm for “climbing up” the circles or spheres. This process, for the circles 
at least, is dual to the Euclidean algorithm: if [0; ai, 02,..., a n \ represents a 
rational number q := a/b, then the Gauss map of it (which is tantamount 
to running the Euclidean algorithm one step) satisfies {1 /q} = [0; 02, ..., a„] 
whereas the parents of C{q) are C(q') and C(q”) where q' = [0; a \,..., a„_ 1 ] and 
q" = [0; ai,..., a n — 1]. Hence, the Euclidean property for Z is equivalent to 
the property that every Ford circle of radius less than 1/2 touches a larger one. 
This is also a consequence of having a recursive geometric procedure (that adds 
smaller circles) which generates all of the Ford circles. 

It is reasonable then to consider the tetrahedral and octahedral geometric 
recursions of Sections 4 and 5 as ways to “see” that Z[i\ and Z[w] are Euclidean 
domains. 

Problem 7.3. It is well known that there are exactly five imaginary quadratic 
fields whose rings of integers are Euclidean. Are there then only five polyhe- 
dra that define a recursive geometric procedure that gives all the Ford spheres 
parameterized by these fields? 

A known theorem states that the integers of an imaginary quadratic field 
form a UFD if and only if it satisfies a multistage Euclidean algorithm (Propo¬ 
sition 3.2 of m)- A celebrated theorem of Heegner, Stark, and Baker (all 
independent) showed that there are precisely 19 such fields. 

Problem 7.4. What polyhedra form contact graphs for spheres and whose it¬ 
erates form all the Ford spheres for the 19 UFD’s? 

Work in this direction seems to have been done by Yasaki ,32] . 

Given any circle packing in C, normal spheres can be attached to each point 
of tangency of these circles (assign curvature to the sphere equal to the sum of 
curvatures of the two circles; see [20j ) so as to form an array of non-overlapping 
normal spheres. By the Koebe-Andreev-Thurston theorem, every planar graph 
can be realized as the contact graph of some circle packing and thus give rise to 
the contact graph of normal spheres. 
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Problem 7.5. Every finite contact graph of spheres generates a recursive geo¬ 
metric process leading to an infinite family of spheres whose points of tangency 
with C form a set with some type of Euclidean property; can this idea be devel¬ 
oped and/or can such a set of be algebraically defined? 

Although we expanded V to V„ in Sections 4,5, and 6, for imaginary quadratic 
fields with class number 1, it seems possible to extend further. 

Problem 7.6. Investigate V a for a, say , a cubic algebraic number. 

Problem 7.7. Solutions of equation (11) parameterized Ford circles V a ; do 
more general equations parameterize other interesting families of spheres ? 

Problem 7.8. A conjecture by Gauss states that there are infinitely many real 
quadratic fields of class number 1 (i.e., ring of integers is a UFD). Are there 
definitions of Ford circles that are relevant to this conjecture? 

For example, Q(-\/2) is a UFD. Given a, b,c,d £ Z[\/2], consider the recursive 
geometric procedure for normal circles: 

{C a ,b, C'c.d} '-t {C a ,b, C ay / 2 +c by / 2 +d , C a+c ^/2 } b+dV2^ ^c,d} 

The points where these circles intersect R take on values that form a strict 
subset of Z[-\/2j. These circles form a subset of a circle packing investigated by 
Guettler and Mallows [8]. 

At the end of Section 6, it was seen Ford spheres parameterized by Q(er) can 
also be parameterized by a group {(x, y, z) £ Q : x ® y ® z = ej where © is the 
“secant addition” m associated with the minimum polynomial for a. 

Problem 7.9. A group of rationals (possibly with e = oo) can be based on secant 
addition for any quadratic polynomial. What families of spheres are generated 
in this way? Is there a way to understand the group structure geometrically 
(i.e., in terms of the spheres)? Considering cubic polynomials in this context, 
© is still defined though not associative. Is there a way to define spheres in this 
case? 
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